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Using the variational principle, we compute mass spectra and decay constants of ground state pseudoscalar 
and vector mesons in the light-front quark model (LFQM) with the QCD-motivated effective Hamiltonian in¬ 
cluding the hyperfine interaction. By smearing out the Dirac delta function in the hyperfine interaction, we avoid 
the issue of negative infinity in applying the variational principle to the computation of meson mass spectra and 
provide analytic expressions for the meson mass spectra. Our analysis with the smeared hyperfine interaction 
indicates that the interaction for the heavy meson sector including the bottom and charm quarks gets more point¬ 
like. We also consider the flavor mixing effect in our analysis and determine the mixing angles from the mass 
spectra of (co,tj)) and ( 77 , 77 '). Our variational analysis with the trial wave function including the two lowest 
order harmonic oscillator basis functions appears to improve the agreement with the data of meson decay con¬ 
stants and the heavy meson mass spectra over the previous computation handling the hyperfine interaction as 
perturbation. 


I. INTRODUCTION 

Effective degrees of freedom to describe a strongly inter¬ 
acting system of hadrons have been one of the key issues in 
understanding the non-perturbative nature of QCD in the low 
energy regime. Within an impressive array of effective theo¬ 
ries available nowadays, the constituent quark model has been 
quite useful in providing a good physical picture of hadrons 
just like the atomic model for the system of atoms. Absorbing 
the complicated effect of quark, antiquark and gluon interac¬ 
tions into the effective constituent degrees of freedom, one 
may make the problem more tractable yet still keep some key 
features of the underlying QCD to provide useful predictions 
[1], The effective potentials used in constituent quark models 
are typically described by the flux tube configurations gener¬ 
ated by the gluon fields as well as the effective “one-gluon- 
exchange” calculation in QCD [2, 3], In the QCD-motivated 
effective Hamiltonian, a proper way of dealing with the rel¬ 
ativistic effects in the hadron system is quite essential due to 
the nature of strong interactions. In particular, proper care 
and handling of relativistic effects has been emphasized in 
describing the hadrons made of it, d, and v quarks and anti¬ 
quarks. 

As a proper way of handling relativistic effects, the light- 
front quark model (LFQM) [4—8] appears to be one of the 
most efficient and effective tools in hadron physics as it 
takes advantage of the distinguished features of the light- 
front dynamics (LFD) [9, 10], In particular, the LFD car¬ 
ries the maximum number (seven) of the kinetic (or interac¬ 
tion independent) generators and thus the less effort in dy¬ 
namics is necessary in order to get the QCD solutions that 
reflect the full Poincare symmetries. Moreover, the ratio¬ 
nal energy-momentum dispersion relation of LFD, namely 
p~ = (p^ +m 2 )/p + , yields the sign correlation between the 
light-front (LF) energy p (= p° - p 1 ’) and the LF longitudi¬ 
nal momentum p + (= p° + p') and leads to the suppression 
of quantum fluctuations of the vacuum, sweeping the compli¬ 
cated vacuum fluctuations into the zero-modes in the limit of 
p + —> 0 [11-13]. This simplification is a remarkable advan¬ 


tage in LFD and facilitates the partonic interpretation of the 
amplitudes. Based on the advantages of the LFD, the LFQM 
has been developed [14] and subsequently applied for vari¬ 
ous meson phenomenologies such as the mass spectra of both 
heavy and light mesons [15], the decay constants, distribu¬ 
tion amplitudes, form factors and generalized parton distribu¬ 
tions [10, 14-23], 

Despite these successes in reproducing the general features 
of the data, however, it has proved very difficult to obtain di¬ 
rect connection between the LFQM and QCD. Typically, rig¬ 
orous derivations of the connection between the effective con¬ 
stituent degrees of freedom and the fundamental QCD quark, 
antiquark and gluon degrees of freedom have been explored 
by solving momentum-dependent mass gap equations as dis¬ 
cussed in many-body Hamiltonian approach [24], Dyson- 
Schwinger approach [25], etc. Although one has not yet 
explored solving the momentum-dependent mass gap equa¬ 
tion in LFD, there has been some attempt to derive an ef¬ 
fective LF Hamiltonian starting from QCD using the discrete 
light-cone quantization (DLCQ) and solve the correspond¬ 
ing equation of motion approximately for the quark and anti¬ 
quark bound-states to provide semianalytical expressions for 
the masses of pseudoscalar and vector mesons [26], The at¬ 
tempt to link between QCD and LFQM is also supported 
by our recent analyses of quark-antiquark distribution am¬ 
plitudes for pseudoscalar and vector mesons in LFQM [27], 
where we presented a self-consistent covariant description of 
twist 2 and twist 3 quark-antiquark distribution amplitudes 
for pseudoscalar and vector mesons in LFQM to discuss the 
link between the chiral symmetry of QCD and the LFQM. 
Our results for the pseudoscalar and vector mesons [27] ef¬ 
fectively indicated that the constituent quark and antiquark in 
the LFQM could be considered as the dressed constituents in¬ 
cluding the zero-mode quantum fluctuations from the QCD 
vacuum. Moreover, the light-front holography based on the 
5-dimensional anti-de Sitter (AdS) spacetime and the confor¬ 
mal symmetry has given insight into the nature of the effec¬ 
tive confinement potential and the resulting light front wave- 
functions for both light and heavy mesons [28]. As we have 
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shown in Ref. [29], our LFQM analysis of the pion form fac¬ 
tor provided compatible results both in spacelike and time¬ 
like regions with the holographic approach to LF QCD [30]. 
These developments motivate our present work for the more- 
in-depth analysis of the mass spectra and decay constants for 
the ground state pseudoscalar and vector mesons in LFQM. 

In LFQM, the LF wave function is independent of all ref¬ 
erence frames related by the front-form boosts because the 
longitudinal boost operator as well as the LF transverse boost 
operators are all kinematical. This is clearly an advantageous 
feature unique to LFQM, which makes the calculation of ob¬ 
servables such as mass spectra, decay constants, form factors, 
etc. much more effective. Computing the meson mass spec¬ 
tra, however, we have previously [14, 15] treated the hyper- 
fine interaction as a perturbation rather than including it in 
the variation procedure to avoid the negative infinity from the 
Dirac delta function contained in the hyperfine interaction. In 
the present work, we smear out the Dirac delta function by a 
Gaussian distribution and resolve the infinity problem when 
variational principle is applied to the hyperfine interaction. 
We obtain optimal model parameters in our variational anal¬ 
ysis including the hyperfine interaction and examine if it im¬ 
proves phenomenologically our numerical results compared 
to the ones obtained by the perturbative treatment of the hy¬ 
perfine interaction. For our trial wave function, we also take a 
larger harmonic oscillator (HO) basis to see if it provides any 
phenomenological improvement in our predictions of mass 
spectra and decay constants for ground state pseudoscalar and 
vector mesons. 

The paper is organized as follows: In Sec. II, we describe 
our QCD-motivated effective Hamiltonian with the smeared- 
out hyperfine interaction. Using the mixture of the two lowest 
order HO states as our trial wave function of the variational 
principle, we find the analytic formula of the mass eigenval¬ 
ues for the ground state pseudoscalar and vector mesons. The 
optimum values of model parameters are also presented in 
this section. In Sec. Ill, we present our numerical results of 
the mass spectra obtained by taking a larger HO basis in the 
trial wave function and compare them with the experimen¬ 
tal data as well as our previous calculations [14, 15]. To test 
our trial wave function with the parameters obtained from the 
variational principle, we also calculate the meson decay con¬ 
stants and compare them with the experimental data as well as 
other available theoretical predictions. Summary and conclu¬ 
sion follow in Sec. IV. The detailed procedure of fixing our 
parameters through variational principle is presented in Ap¬ 
pendix A. 


II. MODEL DESCRIPTION 

As mentioned in the introduction, there has been an attempt 
to derive an effective LF Hamiltonian starting from QCD us¬ 
ing DLCQ [26], Transforming the LFD variables to the ordi¬ 
nary variables in the instant form dynamics (IFD), one may 
see the equivalence between the resulting effective LF Hamil¬ 
tonian for the quark and antiquark bound-states and the usual 
relativistic constituent quark model Hamiltonian for mesons 


typically given in the rest frame of the meson, i.e. the cen¬ 
ter of mass (C.M.) frame for the constituent quark and anti¬ 
quark system. It may be more intuitive to express the effective 
LF Hamiltonian describing the relativistic constituent quark 
model system for mesons in terms of the ordinary IFD vari¬ 
ables. Effectively, the meson system at rest is then described 
as an interacting bound system of effectively dressed valence 
quark and antiquark typically given by the following QCD- 
motivated effective Hamiltonian in the quark and antiquark 
C.M. frame [14, 15]: 

Hc.m. = \J m 2 q + k 2 + \J w? +P + V, (1) 

where k = (k ; ,k z ) is the relativistic three-momentum of the 
constituent quarks and V is the effective potential between 
quark and antiquark in the rest frame of the meson. The ef¬ 
fective potential V is typically given by the linear confining 
potential V con f plus the effective one-gluon-exchange poten¬ 
tial Vo ge . For 5-wave pseudoscalar and vector mesons, the ef¬ 
fective one-gluon-exchange potential reduces to the coulomb 
potential V cou i plus the hyperfine interaction Vh yp . Thus, one 
may summarize V as 

V = V con f + Voge 

coul hyp 

_ I 7 4OJ s , 2 Sq ■ Sq 2 t l 

— a + br --— + —-V Vcoui, (2) 

v - V -' 3 r 3 ItlqUlq 


where a s is the strong interaction coupling constant 1 , (S q ■ 
S q) = 1/4 (—3/4) for the vector (pseudoscalar) meson and 
V 2 V C oui = (167ra v /3)(5 :, (r). Reduction of the LF Hamilto¬ 
nian in QCD to a similar form of the effective Hamiltonian 
in the C.M. frame of the quark and antiquark system given by 
Eqs. (1) and (2) was discussed in Ref. [26]. For the hyper¬ 
fine interaction Vh yp , one may consider the relativization such 
as Vhyp y \JHlqlTlq/EqEqVhyp -\JHlqlllqjEqEq [31, 32]. Such 
relativization may be important for the <5 3 (r)-type potential 
without any smearing in computing particularly the light me¬ 
son sector. Since we apply the variational principle even for 
the hyperfine interaction in this work smearing out the Dirac 
delta function to resolve the infinity problem, we naturally in¬ 
troduce a smearing parameter which may effectively compen¬ 
sate the factor due to the relativization. With this treatment, 
we are able to provide explicit analytic expressions for the 
meson mass spectra (see Eq.(9)). 

While the effective bound-state mass square is given 
by Mqq = (Pq m ) 2 in the C.M. frame of the constituent quark 
and antiquark system, the energy-momentum dispersion re¬ 
lation in LFD is given by = P + P — P 2 , where the 
four-momentum of the bound system is denoted by P^ = 
(P + ,P~, Px) = (P° +P 3 ,P° — P 3 ,Px). From this, one may 


1 Although one may consider a running coupling constant, we take a s as one 

of the variation parameters in this work. 
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consider the LFQM mass square operator P P — P 2 (that 
provides the eigenvalues P + P — P 2 ) as the square of the 
effective Hamiltonian given by Eq. (1), i.e. Hq M - Since 
the eigenvalues and the expectation values are same for the 
eigenstates, we compute the expectation value (Hc.m.) using 
the variation principle. Alternatively, one may consider com¬ 
puting the expectation value (//<?, M ) in view of the LFQM 
mass square operator P + P~ I* 2 being (Hq M ). Although 
((A//c.m.) 2 ) = 0 in principle for the eigenstates, it may be in¬ 
teresting to examine numerically how small the corresponding 
deviation ((A//c.m.) 2 ) = (^cm ) “ (Hc.m.) 2 is. More future 
works complementary to our present computation of (Hc.m.) 
can be suggested in variational analysis. In this work, we ex¬ 
amine the x 2 values of our computational results in compari¬ 
son with experimental data to get optimal parameter values in 
the (Hc.m.) computation. This will provide useful ground in¬ 
formation for any alternative and/or further works beyond the 
present analysis. 

As discussed earlier, the longitudinal boost operator as well 
as the LF transverse boost operators are all kinematical and 
thus the LF wave function does not depend on the external 
momentum, i.e. P + and I’ . In effect, the determination of 
the LF wave function in the meson rest frame such as P + = 
Mqq and P i = 0 won’t hinder its use for any other values of 
P + and Pj_. This provides the applicability of LFQM for the 
computation of observables beyond the meson mass spectra. 

The wave function is thus represented by the Lorentz in¬ 
variant internal variables x,- = pf /P . k l( = px, ■ x,P and 
helicity A,, where /)( L is the momenta of constituent quarks. 
Explicitly, the LF wave function of the ground state mesons is 
given by 


^obfok-iiA) =^X(*i.k±i)« , (*f»k±j)» (3) 

where <t> is the radial wave function and is the 

AqAq 

interaction-independent spin-orbit wave function. The spin- 
orbit wave functions for pseudoscalar and vector mesons are 
given by [14, 33] 


\f2jMl - ( lllq-niq ) 2 


— 

XqX,q 


u \(Pq) t(Jz) Mq+mJ+lq V A?(F<7) 

\/2 JmI - (itlq - Itlq) 2 


where £^(J Z ) is the polarization vector of the vector meson 
and the boost invariant meson mass squared Mq obtained from 
the free energies of the constituents is given by 


M 2 = 



k l+'»g- 

1 —x 


(5) 


The spin-orbit wave functions satisfy the relation 
JA Ki \ = 1 f° r both pseudoscalar and vector 

mesons. 99 99 

To use a variational principle, we take our trial wave func¬ 


tion as an expansion of the true wave function in the HO ba¬ 
sis. We use the same trial wave function expanded with the 
two lowest order HO wave functions <J> = Y.n=i c ntynS for both 
pseudoscalar and vector mesons, where 


0ts(*i,kxi) 


47Z 3 / 4 jdk z - 



( 6 ) 



and /j is the variational parameter. We should note here that 
our LF wave functions (j)„s are dependent on Ml and thus can¬ 
not be factorized into a function of kx; multiplied by another 
function of x,-. In particular, k 1 in Eqs. ( 6 ) and (7) is given 
by k 2 = k 2 + kp where k z = (x — 1 /2)Mq + (m? — m 2 )/ 2 Mo. 

For instance, e - ^ 2 / 2 ^ 2 = e m2 / 2 P~ e ~ M o/%P~ i n the case of equal 
quark and antiquark mass m q = >riq = m. The variable trans¬ 
formation (x, kx) — ► k = (k 1 ,k z ) requires the Jacobian factor 
given by dk z /dx =Mq[\ — (nr q — w?) 2 /Mq]/4x( 1 — x) as one 
can see from Eqs. ( 6 ) and (7). The normalization of the wave 
function is thus given by 



With <f> = Xh-i < 7 ;b„.y, we evaluate the expectation value 
of the Hamiltonian in Eq. (1), i.e. ( < t , |//c.M.| c E > ) which de¬ 
pends on the variational parameter p. According to the vari¬ 
ational principle, we can set the upper limit of the ground 
state’s energy by calculating the expectation value of the sys¬ 
tem’s Hamiltonian with a trial wave function. In our pre¬ 
vious calculations [14, 15], which we call “CJ model”, we 
first evaluate the expectation value of the central Hamilto¬ 
nian T + Vconf + Vcoui with the trial function <j>ig, where T 
is the kinetic energy part of the Hamiltonian. Once the 
model parameters are fixed by minimizing the expectation 
value (<j>is\(T +V C onf + Vc OU i) \(j>is), then the mass eigenvalue 
of each meson is obtained as M q q = ((j>is\Hc.M. |$is). The hy- 
perfine interaction Vh yp in CJ model, which contains a Dirac 
delta function, was treated as perturbation to the Hamilto¬ 
nian and was left out in the variational process that optimizes 
the model parameters. The main reason for doing this was 
to avoid the negative infinity generated by the delta function 
as was pointed out in [31], Specifically, (0is|Vhyp|0is) for 
pseudoscalar mesons decreases faster than other terms that in¬ 
crease as /j increases and the expectation value of the Hamil¬ 
tonian is unbounded from below. 

The singular nature of the hyperfine interaction and its reg¬ 
ularization is a standard topic in atomic physics and the atomic 
analysis has been carried out to extraordinary precision [34]. 
In particular, a Bethe-Salpeter based bound-state formalism 
was applied to the calculation of recoil contributions of order 
mdP to hyperfine splitting in ground-state positronium [35]. 
Instead of dropping the relative energy dependence in favor 
of equations with a simpler kinematical structure but a more 
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complicated effective kernel, the Barbieri-Remiddi formal¬ 
ism [36] was discussed as an effective way to handle signif¬ 
icant complications concerning the Bethe logarithm [37]. As 
discussed in Ref. [35], the 8 function of the relative energy 
po is replaced by a smearing function of po in the Barbieri- 
Remiddi formalism [36]. In LFD, the equal LF time x + (= 
x° + x 3 * ) correlates the ordinary time x° and space x 3 so that 
the idea of smearing p° in the Barbieri-Remiddi formalism 
may be extended to smear the (5 3 (r) function in hyperfine in¬ 
teraction discussed in the present work. In this respect, our 
regularization procedure discussed below would also be valid 
and is compatible for the hyperfine splitting in atoms. Ana¬ 
lytic treatment of positronium spin splittings was presented in 
LF QED [38] and more recent DLCQ application to the anal¬ 
ysis of H + ll~ bound state spectrum can be found in Ref. [39]. 

To avoid the negative infinity, we thus use a Gaussian 
smearing function to weaken the singularity of 5 3 (r) in hy¬ 
perfine interaction, viz. [31, 32], (5 3 (r) —► (o 3 /jl 3 / 2 )e~ a ~ T ~. 
Once the delta function is smeared out like this, a true min¬ 
imum for the mass occurs at a finite value of p. The ana¬ 
lytic formulae of mass eigenvalues for our modified Hamilto¬ 
nian with the smeared-out hyperfine interaction, i.e. M qq = 
( < J>|f/c.M. | < t > ), are found as follows: 2 




— Q. - h 


b 



JL 

y/n 


£ | \[k ^\/ 6 ciC 2 — 3c 2 ^ U 





+ l c 2 (3 -Zi)K 2 (j) 

+ jUe^ (2c?z 2 —3q —6%/6 c 1 c 2 +9) K x (^)j 

- |5+c 2 + 6\/273ciC2 


4p 2 a 3 (S q -S q ) 
(/3 2 + <7 2 ) 1/l: m q m q 


(2\/6ciC2 + 3-c{ S j 


+2p 2 ^ 2 c 2 + v / 6 ciC 2 ^j 



(9) 


where Zi = »t 2 /j3 2 and K\ is the modified Bessel function of 
the second kind and U(a,b,z) is Tricomi’s (confluent hyper¬ 
geometric) function. We should note that the mass formula for 
the delta-function hyperfine interaction corresponds to Eq. (9) 
in the limit of a —> °°. We then apply the variational princi¬ 
ple, i.e. dM qq /dp = 0, to find the optimal model parameters 
in order to get a best fit for the mass spectra of ground state 
pseudoscalar and vector mesons (a more detailed description 
of this procedure can be found in A). 

Our optimized potential parameters are obtained as {a = 
-0.6699 GeV,/? = 0.18 GeV 2 , a s = 0.4829}. For the best fit 


2 Although the true minimum occurs with the smeared-out hyperfine interac¬ 

tion even for (j)\s case, we found that the phenomenological results do not 

show any significant improvement compared to CJ model. 


of the ground state mass spectra, we obtain ci = + \/0J and 
C 2 = +V0.3. We should note that our potential parameters 
are quite comparable with the ones suggested by Scora and Is- 
gur [40], where they obtained a = —0.81 GeV, b = 0.1 8 GeV 2 , 
and a s = 0.3 ~ 0.6. For a comparison, the coupling constant 
we found in our previous CJ model [14, 15] was a s = 0.31. 

While we use the common potential parameters (a,b,a s ) 
for all the mesons, it was shown in [32, 41] that if a smearing 
procedure for the 5 3 (r) function is used, then a large Gaus¬ 
sian parameter cr is obtained for the heavy quark sector. In our 
updated potential model using the smeared hyperfine interac¬ 
tion (<7 3 /7T 3 // 2 )e“ crr , we also confirm the same observation 
as in [32] for the heavy meson sector including ( b,c ) quarks. 
Thus, we differentiate the smearing parameter cr for the heavy 
(b,c) sectors such as (cc,bc,bb) from the other (qq) sectors by 
introducing multiplicative factor in front of cr, i.e. cr —r Act 
with A > 1 , while other potential parameters (a,b,a s ) remain 
the same for all (qq) sectors. This differentiation is to ac¬ 
commodate the hyperfine splittings for the heavy (h.c) quark 
sectors as we will show in the next section. Our new updated 
results with A differentiation using the common potential pa¬ 
rameters (a,/?, 0^,(7) for all meson sectors show definite im¬ 
provement in the y 1 fit of the experimental data for meson 
masses. 

Our optimal constituent quark masses and the smearing pa¬ 
rameters cr are listed in Table I. Since we included the hyper¬ 
fine interaction with smearing function entirely in our varia¬ 
tional process, we now obtain the two different sets of P val¬ 
ues, one for pseudoscalar and the other for vector mesons, 
respectively. The optimal Gaussian parameters p qq for pseu¬ 
doscalar and vector mesons are also listed in Table II. We 
should note that the values of the multiplicative factor A to 
get the best fits for the mass eigenvalues are obtained as 
A = (2,2.3,3) for (cc,bc,bb) sectors. As a sensitivity check, 
however, we present the numerical results with the follow¬ 
ing theoretical error bars A = ( 2 +}, 2 . 3 l}, 3 + 2 ) for (cr, bc.bb) 
sectors, respectively. Although one may fine-tune more to im¬ 
prove the hyperfine splittings for the heavy-light sectors by 
using different set of A parameters, we set A = 1 for any other 
qq sectors except (cc, be. bb) sectors in this work for simplic¬ 
ity. 


TABLE I. Constituent quark masses [GeV] and the smearing parame¬ 
ter a [GeV] obtained by the variational principle for the Hamiltonian 
with a smeared-out hyperfine interaction. Here q = u and d. 


m q 

m s 

m c 

m b 

a 

0.205 

0.380 

1.75 

5.15 

0.423 


We also determine the mixing angles from the mass spec¬ 
tra of (o,0) and ( 77 , 77 '). Identifying (&,&') = {(jf.co) and 
( 77 , 77 ') for vector and pseudoscalar nonets, the flavor assign¬ 
ment of & and mesons in the quark-flavor basis ml = 
(uu + dd)/\/l and ss is given by [42-44] 

( & \ ( cosa —sina \ ( nh \ TT , . / nn 

(jr ) = ( sina cosa ) ( ss }= U <“>(**)W 

where a is the mixing angle in the quark-flavor basis. For 
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TABLE II. The Gaussian parameter j8 [GeV] for ground state pseudoscalar (J pc = 0 l ") and vector (1 ) mesons obtained by the variational 

principle, q = u and d. We should note that A = (2^],2.3+[, 3f 2 ) are used to get (pec, PbcPbb) values and A = 1 is used to get the rest of j8 qq 
values. 


J PC 


Ais 

j3 ss 

/)qc 

Pcs 

Pcc 

Pqb 

Pbs 

Pbc 

Pbb 

0-+ 

0.4465 

0.3759 

0.3445 

0.3801 

0.3859 

0 5270 +U U291 

u.jz/u_ 00235 

0.4226 

0.4412 

0 6646 +U U219 

u.oo4fo_ 00174 

0 9906 +U U42U 

u.^uo_ 00223 

1 

0.2346 

0.2598 

0.2820 

0.3445 

0.3667 

0 4Ql4 u ut,b/! 

U.H-71H- q0 0072 

0.4057 

0.4321 

n 0.0058 

0.6365 +aoo56 

n 0.0122 

0.9603 0(Km 


the rj — 77 ' mixing, the SU(3) mixing angle 0 in the the fla¬ 
vor SU(3) octet-singlet basis (rjg.tji) can also be used and 
the relation between the mixing angles is given by 0 = a — 
arctan\/2 ~ a — 54.7°[45], Taking into account SU(3) sym¬ 
metry breaking and using the parametrization for the (mass) 2 
matrix suggested by Scadron [46], we obtain [14] 


tan 2 a = 




( 11 ) 


which is the model-independent equation for any qq meson 
nonets. The details of obtaining meson mixing angles using 
quark-annihilation diagrams are summarized in [14], where 
the mixing angle 8 = a — 90° is used in the quark-flavor 
basis. In order to predict the co — (j) and /] T]' mixing 

angles, we use the experimental values of Mjr = (M^,Mjj) 
and M^ri = (M m ,M q i) as well as the masses of M^~ [M^-\ = 
780 (901) MeV and M p n [Af£] = 140 (726) MeV obtained 
from (<t>|// s .s|<t>) for both vector (V) and pseudoscalar (P) 
mesons, respectively. Our prediction for co — (j) mixing an¬ 
gle is = 84.8°, which is about 5.2° deviated from the 
ideal mixing = 90°. Our prediction for 1 ] — 1 ]' mixing 
angle is a r| _^ = 36.3°, which is in agreement with the range 
34.7° to 44.7° of phenomenological values [42, 45]. 

Our updated model with the smeared hyperfine interaction 
appears to improve the result of mass spectrum, which is pre¬ 
sented in the next section. This may suggest that when using 
constituent quark models, the contact interactions has to be 
smeared out in general. In fact, we think this smeared interac¬ 
tion is more consistent with the physical picture for a system 
of the effective constituent quarks which are not point-like. 

For practical application of our model, we also compute the 
decay constants for the ground state pseudoscalar and vector 
mesons. The decay constants are typically defined by 


{0\qy^y 5 q\P)=ifpP^ 

(0\qfq\V(P,h))=fvM v e fl (h), 


for pseudoscalar and vector mesons, respectively. For the /] 
and 17 ' case, one may also define decay constants through ma¬ 
trix elements of octet and singlet axial-vector currents. How¬ 
ever, as discussed in [42, 43], they cannot be expressed as 
U (0)diag[/g,/i] due to the £7(1)^ anomaly. Thus, the follow¬ 
ing two mixing angle parametrization is adopted [42, 43] 

4 = h cos 08, 4 = ~fi sin 0i, 

4, =/g sin 0 8 , 4,=/! cos 0’. (13) 

The parameters appearing in Eq. (13) are related to the basis 


parameters a,f q = f n n and f s = f sS , characterizing the quark- 
flavor mixing scheme as follows [42]: 


, 2 4+ 2 4 

h = —?— 


, 2 2/ 2 + / 2 


Gg = a — arctan 


fi = 


01 = a — arctan 


y/2 fs ) 

U ) 


tm. (. 4 ) 


Using the plus component (/i = +) of the currents, one can 
calculate the decay constants. The explicit formulae of pseu¬ 
doscalar and vector meson decay constants in quark-flavor ba¬ 
sis are given by [14, 33] 


fp 

fv 



d 2 kx ^ 

’ 

d 2 ki <4>(x,k I ) [ 2k 2 

-r- , = - 4 , 

8?r ^2 + k 2 I- DLF J 


(15) 


where = (1 — x)m 9 + xm q and D\ \ = Mo + m q + m q . 


III. RESULTS AND DISCUSSION 


In Fig. 1, we show the masses (and hyperfine splittings) 
and the corresponding decay constants of heavy quarkonia de¬ 
pending on the variation of the multiplicative factor A. For 
(cc) (Fig. 1(a) and 1(b)) and ( bb) (Fig. 1(c) and 1(d)) , we 
plot the curves corresponding three different A values, i.e. 
Act = (1,2,3)a and Act = (1,3,5)cr, respectively. The results 
indicate that the interaction between heavier quarks gets more 
point-like as the larger A values are favored in comparison 
with data. In general, as one can see from Fig. 1, the hyper¬ 
fine splittings for both charmoninum (Fig. 1(a)) and bottomo- 
nium (Fig. 1(c)) states increases as A increases while other 
potential parameters remain the same. On the other hand, the 
decay constants of vector mesons decrease while the corre¬ 
sponding decay constants of pseudoscalar mesons increase as 
A increases (see Fig. 1(b) and 1(d)). One may increase A value 
even further to have better hyperfine splittings compared to the 
data. However, one may not increase A value arbitrarily to ac¬ 
commodate the empirical constraint fv > fp. Similarly, we 
could improve the hyperfine splitting for (be) sector by using 
A =2.3. 

We show in Fig. 2 our prediction of the meson mass spectra 
obtained from the variational principle to the effective Hamil¬ 
tonian with the smeared-out hyperfine interaction using the 
trial function 4> = £ 2 c,0,v and compare them with the exper- 
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FIG. 1. (color online). The masses (and hyperfine splittings) and the corresponding decay constants of heavy quarkonia depending on the 
variation of the multiplicative factor A, i.e. (cc) (Fig. 1(a) and 1(b)) with two different A a = (1,2)<r and ( bb ) (Fig. 1(c) and 1(d)) with two 
different A <J = (1,3)<7 values, respectively. 
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FIG. 2. (color online). Fit of the ground state meson masses [MeV] 
with the parameters given in Table II and I compared with the fit 
from our previous calculations using CJ model [15] as well as the 
experimental values. The (jr,p) masses are our input data. The 
(r/,r/' ,£0,0) masses are also used as input to find the (rj — r )') and 
(to — <j>) mixing angles. The theoretical error bars for (cc,bc,bb) 
sectors are due to the usage of A = (2+[ ,2.3^[, 3^,) values, respec¬ 
tively. 


imental data [45] as well as the results obtained from the CJ 
model with the linear confining potential [14]. We should note 
that the (n,p) masses are used as inputs. The (rj,r) , ,®,0) 
masses are also used as inputs to find the (rj — t] 1 ) and (co — p) 
mixing angles. The theoretical error bars for (cc.bc.hh) sec¬ 
tors are due to the usage of A = values, re¬ 

spectively. As one can see, our trial wave function <t> includ¬ 
ing more HO basis generates overall better results than our 
CJ model. This can be seen from our x 2 = 0.008 compared to 
X" = 0.012 obtained from the CJ model [15]. Except the mass 
of K, our predictions for the masses of 15-state pseudoscalar 
and vector mesons are within 4% error. Especially, our effec¬ 
tive Hamiltonian with the smeared hyperfine interaction using 
<f> clearly improves the predictions of heavy-light and heavy 
quarkonia systems such as (r\ c ,J /\jf, B c . f]/,, T) compared to 
the CJ model adopting the contact hyperfine interaction. Al¬ 
though the experimental data for B* is not yet available, our 
predictions of B i.e. 6330 1 ) MeV, are quite comparable 
with the lattice prediction 6331(9) MeV [47] as well as other 
quark model predictions such as 6340 MeV [32] and 6345.8 
MeV [48], 

In Table III, we list our predictions for the decay constants 


TABLE III. Decay Constants for light mesons (in unit of MeV) ob¬ 
tained from our updated LFQM. 


Model 

fn 

U 

.h 

fK* 

This work 

130 

205 

161 

224 

CJ [16] 

130 

246 

161 

256 

Exp. [45] 

130.4(2) 

208 (fl ),216(5)W 

156.1(8) 

217(7) 


W Exp. value for T(t -4 pv T ) 
( h' 1 Exp. value for p° —> e + e~. 


TABLE IV. Decay constants in the singlet-octet basis and the mixing 


angle in the quark-flavor basis. 

Reference 

fs/ fn 

08 

fl/fn 

01 

a 

This work 

1.30 

-27.3° 

1.16 

-8.6° 

36.3° 

[42] 

1.26 

-21.2° 

1.17 

-9.2° 

39.3° 

[44] 

1.28 

-20.5° 

1.25 

-4° 

— 

[49] 

1.51 

-23.8° 

1.29 

-2.4° 

40.7° 

[50] 

1.27 

-19.5° 

1.17 

-5.5° 

42.1° 


of light mesons {n.K.p,K*) obtained by using the mixed 
wave function 4> of 15 and 25 HO states and compare them 
with the results from the CJ model [16] and the experimen¬ 
tal data [45]. As one can see, our updated model calculation 
including the hyperfine interaction in the variation procedure 
clearly improves the results over the CJ model. 

For the decay constant of <j> meson, our prediction for the 
ideal mixing angle = 90°) is given by fy = = 245.1 

MeV. However, we obtain f, = fj s = 226 MeV using our pre¬ 
dicted mixing angle a a ^if, = 84.8°. Comparing to the exper¬ 
imental value /“ p ' = 233 MeV [45] (extracted from the par¬ 
tial width of (j) —> e + e~ decay), our prediction for fo prefers 
a rather small CO — (j) mixing angle such as cx m -<]> ~ 87.5° than 
the ideal mixing. 

For the decay constants of T] and rj\ our predictions of the 
decay constants f q and f s are given by f q = 130 MeV and 
f s = 184.8 MeV so that f q lf K = 1 and f s /f% = 1.42, where 
the SU(3) breaking effect is manifest in the ratio f q /f s 1. 
Using Eq. (14), we obtain f&/f K = 1-30 and fi/fn = 1.16 
with 08 = —27.3° and 9i = —8.6°, respectively. In Table IV, 
we compare our results for the decay constants in the singlet- 
octet basis and the mixing angle in the quark-flavor basis with 
other theoretical predictions [42, 44, 49, 50]. As one can see, 
our results are consistent with other theoretical model results. 
Since the experimental values are very well known for light 
mesons, this improvement is very encouraging. 

In Table V, we list our predictions for the charmed me¬ 
son decay constants (f D , f D *, f Ds , f D *, f Ac Jj^) together with 
CJ model [23], lattice QCD [51-54], QCD sum rules [55], 
relativistic Bethe-Salpeter (BS) model [56], relativized quark 
model [57], and other relativistic quark model (RQM) [58] 
predictions as well as the available experimental data [45, 59]. 
We extract the experimental value (/y/>p)exp = (407±5) MeV 
from the data r exp (T/’P —>• e+e~) = (5.55 ±0.14) keV [45] 
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and the formula 

r(V -A e+e ) = ~Y^QED e Q^^ (16) 

where eg is the electric charge of the heavy quark in units of e 
(2/3 for c and —1/3 for b). We should note that our results of 
the ratios f D Jfo = 1.11 and fnJfj/y = are quite 

comparable with the available experimental data, ZdJZd = 
1.25 ±0.06 [45] and fj/xy = 0.81 ±0.19 [45, 59], respec¬ 
tively. Our result of the ratios fp>* / fry = 1.13 is also in good 
agreement with other theoretical model calculations such as 
1.16 ±0.02 ±0.06 from the lattice QCD [54] and 1.10 ±0.06 
from the BS model [56]. 

We list our results for the bottomed mesons 
C ) in Table VI, and compare with 
CJ model [23], lattice QCD [51, 52, 60], QCD sum rules 
[55], BS model [56], relativized quark model [57], and 
RQM [58] predictions as well as the available experi¬ 
mental data [45, 61], Note that we extract the exper¬ 
imental value (/x)exp = (689 ± 5) MeV from the data 
r exp (T -A e+e~) = 1.340 ±0.018 keV [45] and Eq. (16) 
with e 2 Q = 1/9 for V = Y. Our results for the ratios 
Zb s /Zb = 113 and Zb* /Zb* = 1.15 are in good agreement 
with the QCD sum rules [55] predictions: ZbJZb = 1.17+q^, 
and Zb* I Zb* = 1.20 ±0.04. Ours are also in good agree¬ 
ment with the lattice results, Zb s /Zb = 1-206(24) [52] 
and Zb*/Zb* = 1.17(4 )^3 [51], Our result for the ratio 
hJh = 0-99lg;S is consistent with the heavy quark 
symmetry fr\ b /fx = 1 [62], One can also see that for heavy 
charmed and bottomed mesons, the trial wave function (j )b 
produces better results when compared with the experimental 
data as well as the lattice results. 

In Table VII, we present our model predictions for the de¬ 
cay constants of fg c and fa- , and compare them with other 
model calculations [15, 63-68]. Our results are comparable 
with other model calculations. 


IV. SUMMARY AND CONCLUSION 

In this work, we updated our LFQM by smearing out the 
Dirac delta function in the hyperfine interaction to avoid the 
issue of negative infinity in applying the variational principle 
to the computation of meson mass spectra, while our previ¬ 
ous modeKCJ model) used the perturbation method to handle 
the delta function in the contact hyperfine interaction. Using 
the mixed wave function <f> of I .S' and 2 S HO states as the 
trial wave function, we calculated both the mass spectra of the 
ground state pseudoscalar and vector mesons and the decay 
constants of the corresponding mesons. The flavor mixing 
effect has also been implemented for the meson systems of 
(®, 0 ) and (77,T]'). 

The variational analysis with <t> seems to improve the agree¬ 
ment with the data of meson decay constants over the results 
of the CJ model. It also appears to provide the better agree¬ 
ment with data in the heavy meson mass spectra. Accommo¬ 
dating the empirical constraint fy > fp, we have shown that 


the mass spectra and the hyperfine splittings for heavy (b,c) 
quark sector get improved by introducing the multiplicative 
factor A in front of the smearing parameter cr, i.e. a —> A <7, 
in the smeared hyperfine interaction (c^ / n^- 2 )e a T . Our re¬ 
sults indicate that the interaction between heavier quarks gets 
more point-like as the larger A values are favored in compari¬ 
son with data. The distinction between the heavy meson sec¬ 
tor and the light meson sector is rather natural in our LFQM 
analyses. To get more definite conclusion in this respect, fur¬ 
ther analysis of other wave function related observables such 
as various meson elastic and transition form factors may be 
useful. It may be also interesting to analyze radially excited 
meson states using the larger HO basis. 
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Appendix A: Fixation of the model parameters using variational 
principle 


In our model, we assumed SU(2) flavor symmetry and 
have the following parameters that need to be fixed: con¬ 
stituent quark masses (m,,^, m s , m c ,nib), potential parameters 
(a.b ,«/), gaussian parameter /j, and the smearing parameter 
a. For our trial wave function <t> = Y.n=i c ntynS’ we also have 
the mixing factor c„(n =1,2) that we have to adjust. Notice 
that the p values here are not only different for different quark 
combinations, but also different for pseudoscalar and vector 
mesons of the same quark combination. The reason for this 
is that the hyperfine interaction we included in our parame¬ 
terization process gives different contributions to the masses 
of pseudoscalar and vector mesons and thus induces different 
parameterizations under variational principle. 

We now illustrate our procedure for fixing these parameters. 
The variational principle gives us one constraint: 


d_mm = d_M 3l = 

dp dp 


(Al) 


We can use this equation to rewrite the coupling constant a s 
in terms of other parameters and plug it back into Eq. (9) 
and thus eliminate a s . The string tension b is fixed to 
be 0.18 GeV, a well known value from other quark model 
analysis [32, 40, 69], We will leave the quark masses and 
smearing parameter a and the mixing factor c 1 as exter¬ 
nally adjustable variables. We picked a set of values for 
{m u <d)i m s,mc,mb,0,c 1 ) and proceed with the following pro¬ 
cedure to solve for the rest of parameters. 

We are left with 3 more parameters (a,/3^,j3^) for mesons 
of a certain quark combination (qq), where p^, p^ are the 
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TABLE V. Charmed meson decay constants (in unit of MeV) obtained from our updated LFQM. The theoretical error bars for j y ^ come 

from the variation of the smearing parameters (7, i.e. 


Model 

/d 

fir 

fo s 

fo* 

fr\ c 

fj/y 

This work 

CJ [23] 

Lattice [51] 

QCD [52, 53] 
Sum-rules [55] 
BS [56] 

QM [57] 

RQM [58] 

Exp. 

208 

197 

211 ± 3 ± 17 

208 ±7 [52] 

201 +41 

230 ±25 

240 ±20 

234 

206.7±8.9 [45] 

230 

239 

245 ± 20^2 

242^j® 

340 ± 23 

310 

231 

232 

231±12+f 

250 ±7 [52] 
238^23 

248 ±27 

290 ±20 

268 

257.5 ±6.1 [45] 

260 

273 

272 ± 16^2 0 

293+}® 

375 ± 24 

315 

353+^7 

326 

387 ±7 [53] 

292 ±25 

335 ±75 [59] 

361 +7 

360 

418±9 [53] 

459 ±28 

407 ±5 [45] 


TABLE VI. Bottomed meson decay constants (in unit of MeV) obtained from our 
from the variation of the smearing parameters tr, i.e. 

updated LFQM. The theoretical error bars for frj b (T) come 

Model 

/b 

/b* 

//is 

fs* 

frib 

/t 

This work 

CJ [23] 

Lattice [51] 

QCD [52,60] 
Sum-rules [55] 

BS [56] 

QM [57] 

RQM [58] 

Exp. 

181 

171 

179 ± 18+;j 4 

189 ±8 [52] 
207+g 7 

196 ±29 

155 ±15 

189 

229+2^+37 [61] 

188 

185 

196 ±24+2® 

2 10+1 ® 

238 ± 18 

219 

205 

205 

204 ± 16/q 1 
228 ±8 [52] 
242/j l 
216 ±32 
210±20 
218 

216 

220 

229 ± 20/jg 

251+ 14 

272 ± 20 

251 

605t$ 

507 

529 

649 ±31 [60] 

498 ± 20 

689 ±5 [45] 


TABLE VII. Bottom-charmed meson decay constants(in unit of MeV) obtained from 
come from the variation of the smearing parameters O', i.e. /B r (B*)(2.3<r/®). 

our updated LFQM. The theoretical 

error bars for f Bc ^ B , j 

Model This work 

CJ [15] 

[63] [64] 

[65] 

[66] 

[67] 

[68] 

~h c 389/f 

349 

360 433 

500 

460 ± 60 

517 

410± 40 

Tv 391^ 

369 

503 

500 

460 ± 60 

517 

- 


gaussian parameters for pseudoscalar and vector mesons, re¬ 
spectively. Using the masses of K and p as our input values 
for M qq in in Eq. (9), and the condition that our coupling con¬ 
stants a s are the same for all these ground state pseudoscalar 
and vector mesons, we can fix the three model parameters 
(a,Pqq,P qq ) for q = u or d from the following three equations: 


— 0.140, 

(A2a) 

Mp(fiqq-a) = 0.780, 

(A2b) 

a s(fiqq- a ) = CC s {pV-,a). 

(A2c) 


Solving these equations not only gives us the remaining pa¬ 
rameters a.fiandji qq , but also the coupling constant a s 
which we assumed to be the same for all the mesons we con¬ 
sider here. We can then solve for the fi values of all the other 


mesons using the known a s value, by equating the a s expres¬ 
sions for different mesons that we got from Eq. (Al). We thus 
fixed all parameters for the ground state pseudoscalar and vec¬ 
tor mesons we consider here. 

We then assign a different set of values to the externally 
adjustable variables, i.e. (m u ^,m s ,m c ,mb,a,c i), and repeat 
the above procedure until we find a set of values that give best 
fit for the meson mass spectra. 

Through our trial and error type of analysis, we found 
m q = 0.205 GeV. m v = 0.38 GeV,m c = 1.75 GeV./n,, = 
5.15 GeV, a = 0.423 GeV,ci = y/OJ gives best fit. We then 
determine the mixing angles from the mass spectra of (w,0) 
and ( 77 , 77 ') using Eqs.(10) and (11) as we have described in 
Sec. II. In addition, the multiplicative factor A in front of the 
smearing parameter a for the (cc, be. bb) systems were ad¬ 
justed utilizing the hyperfine splittings of (b.c) quark sectors 
as we have discussed in Sec. III. The updated p values with 
this A adjustment are listed in Table II. 
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